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Introduction 

The  long-term  objective  of  this  effort  is  to  develop  a  generic  gene  shuffling-based  technology 
to  rapidly  screen  DNA  libraries  encoding  mammalian  paraoxonase  1  (PONl),  and  thus  to 
identify  variants  that  can  intercept  both  existing  and  emerging  organophosphate-based  chemical 
warfare  nerve  agents  (CWNA).  Enzymes  identified  in  these  screens  should  be  capable  of 
catalytically  neutralizing  the  target  agent  under  physiological  conditions,  thereby  providing  a 
basis  for  development  of  a  new  generation  of  therapeutic  agents  against  CWNA,  in  either  pre-  or 
post- treatment  modality.  The  major  milestone  is  to  integrate  established  components  of 
enhanced  molecular  evolution  techniques  so  as  to  provide  a  means  of  screening  libraries  with 
both  greater  sensitivity  and  higher  throughput.  This  effort  includes  the  detoxification  of  OP- 
based  pesticides 

G-  and  V-type  type  nerve  agents:  tabun  (GA),  sarin  (GB),  soman  (GD),  cyclosarin  (GF)  and 
VX,  are  highly  toxic  organophosphates  (OPs)  that  penetrate  the  body  both  by  inhalation  and  skin 
absorption.  Commercial  OP  pesticides  such  as  parathion,  chlopyrifos,  phorate,  and  malathion, 
are  also  considered  as  threat  chemical  agents,  especially  to  civilian  populations.  In  the  last  6 
years  we  have  used  a  strategy  that  utilizes  enhanced  evolution,  protein  engineering  and  3D 
structural  analysis,  together  with  a  specific  interception  screening  protocol,  to  generate  large 
libraries  of  mammalian  PONl  variants  that  could  be  expressed  in  E.  coli  and  screened  for  their 
ability  to  hydrolyse  OP  nerve  agents.  Structural  analysis  was  then  used  to  direct  the  mutagenesis 
of  relevant  active-site  positions  so  as  to  obtain  highly  effective  bioscavengers.  The  combined 
approach  permitted  the  isolation  of  mutants  with  values  for  the  detoxification  of  in-situ- 

generated  G  agents,  such  as  GD  and  GF,  that  approached  levels  required  for  them  to  be 
considered  as  efficacious  catalytic  bioscavenger  drugs  (>10^  M'mm').  In  vivo  experiments 
conducted  at  ICD  showed  that  several  variants  generated  by  our  group  {e.g.,  4E9,  2D8,  0C9,  and 
VII-D-11)  conferred  full  protection  in  pretreated  guinea  pigs  at  doses  as  low  as  1  mg/kg  protein 
against  IxFDjg  doses  (sc)  of  GA,  GB,  GD,  and  GF.  To  date,  the  best  variants  hydrolyze  GD  and 
GF  at  3.0-5.0x10^  M'mm'  and  GA  and  GB  at  ~7-10-fold  lower  rates.  Key  variability  positions 
at  which  changing  the  residue  produces  improved  hydrolysis  of  G  agents  are  69,  115,  134,  136, 
222,291  and  332. 

Progress  for  V  agents  has  been  slower  than  for  G  agents,  due  both  to  the  lack  of  a  suitable 
surrogate  molecule  that  could  mimic  the  methylphosphonyl(A,A-dialkylamino  ethylthiolate) 
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moiety  of  V-type  agents,  and  to  the  fact  that  under  current  screening  conditions  the  starting 
point,  wild-type  (wt)  rePONl,  does  not  hydrolyze  the  V-type  surrogate,  amiton, 
[(C2H50)2P(0)SCH2CH2N(C2H5)2],  at  a  measurable  rate,  certainly  not  at  a  level  detectable  in 
crude  bacterial  lysates.  However,  variants  in  recent  libraries  (VXR4  and  VXR5),  screened 
directly  with  VX,  displayed  values  of  ~ lx  10^  M" 'min' ',  and  could  be  used  as  lead  mutants 

for  advanced  directed  evolution  of  antidotes  against  both  VX  itself  and  several  VX  analogs,  such 
as  VR. 

Body 

I.  General  Aims 

The  overall  specific  aims  of  the  project  research  at  the  Weizmann  Institute  are: 

1.  Development  of  high-throughput  assays  for  OP  hydrolase  variants  exhibiting  high 
specificity  factors  and  turnover. 

2.  Provision  of  proof-of-concept  for  the  proposed  core  technology  employing  directed 
evolution  of  new  recombinant  PON  and  AChE  variants. 

3.  Isolation  of  interceptors  for  OP-based  commercial  pesticides  and  G-  and  V-type  nerve 
agents,  and  expression  in  soluble  form. 

4.  Large-scale  production  of  selected  enzyme  candidates,  and  their  kinetic,  structural  and 
pharmacological  evaluation 

6.  Establishment  of  “off-the-shelf’  libraries  for  rapid  identification  of  antidotes  against  future 
emerging  threats 

II.  Year  07  specific  milestones  (of  the  renewal  NIH  project): 

The  tasks  to  be  performed  at  the  Weizmann  Institute  to  achieve  the  major  goals  of  the 
Center  for  the  current  year  are  outlined  below: 

#1.  Screening  for  paraoxon-,  parathion-,  malaoxon-,  malathion-,  and  phorate-hydrolyzing 
variants 

#2.  Validation  of  the  hypothesis  concerning  the  structural  motifs  in  PONl  that  are  involved  in 
association  with  HDL,  improving  rePONl ’s  affinity  for  HD,  and  increasing  the  biological  life 
time  of  rePON  variants 
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#3.  Use  of  steered  molecular  dynamics  to  study  the  trajectory  of  interaction  of  VX  with  re  PONl, 
improvement  of  the  modeling  of  VX  into  the  rePONl  scaffold,  and  studies  of  the  mechanism 
of  action  of  PONl  variants  on  V-type  nerve  agents 

#1:  Screening  for  paraoxon-,  parathion-,  malaoxon-,  and  phorate- 
hydrolyzing  variants 

To  obtain  a  lead  for  newly  designed  PON  Is  capable  of  hydrolyzing  effectively  both 
parathion  and  paraoxon,  we  first  tested  several  variants  selected  from  libraries  generated  towards 
G-  and  V-type  nerve  agents.  Evaluation  was  based  on  monitoring  the  release  of  p-nitrophenol. 

Results  obtained  thus  far  have  not  revealed  a  clear  trend  for  parathion.  Screening  with  the 
0X0  metabolite,  paraoxon,  turned  out  to  be  quite  successful  when  libraries  VXG6  and  VXG7 
were  searched  for  mutants  with  the  capacity  to  hydrolysis  it.  Results  are  summarized  in  Table  1. 


Table  1:  Catalytic  activity  towards  paraoxon  of  variants  in  libraries  VXG6  and  VXG7 
(kj,^yKM,  xlO®  M'mm';  figures  in  parentheses  are  SD,  n=2-4). 


Library 

Variant 

kcat/KM?  paraoxon 

VXG6 

637C1 

5.79±0.10 

610H5 

4.44±0.042 

617F7 

3.66±0.39 

6bG12 

1.23±0.06 

617F4 

1.20±0.25 

VXG7 

114-2-Cl 

21.2±2.6 

167-mix-A6 

2.10±0.30 

268-2-E8 

4.69±0.60 

223-1-C3 

1.83±0.12 

167-mix-F5 

1.65±0.18 

As  pointed  out  in  our  previous  (6“'  year)  annual  report,  sub-libraries  of  variants  of  the  6* 
generation  evolved  to  hydrolyze  V  agents  (VXG6),  were  constructed  based  on  the  9  most 
catalytically  efficient  variants  evolved  in  the  5*  round  (VXG5)  of  directed  evolution  for  V-agent 
hydrolysis.  In  addition,  these  were  shuffled  with  oligonucleotides  encoding  for  particular  site- 
specific  mutations.  Specifically,  we  targeted  residues  72,  189,  192,  193,  196  and  292,  which  our 
docking  model  of  VX  with  PONl  implicated  in  substrate  binding,  and  especially  in  leaving- 
group  binding.  In  addition,  we  targeted  residues  69,  71,  73,  74,  115,  134,  222,  291  and  346.  The 
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7th  round  PONl  libraries  of  variants  for  V-type  hydrolysis  (VXG7)  were  generated  by  targeting 
mutagenesis  to  residues  that  flank  both  sides  of  essential  active-site  residues  (E53,  W115,  N168, 
M222,  N224,  D269,  S332),  which  are  involved  in  catalysis  and  with  which  the  catalytic  Ca"^^ 
interacts.  The  underlying  rationale  was  that  changes  in  the  closest  neighbors  of  residues  essential 
for  catalysis  might  help  to  overcome  steric,  geometric  and  electronic  constraints  that  reduced 
catalytic  efficiency  in  the  variants  generated  in  the  earlier  rounds.  These  permutations  were 
examined  vs  the  OP  commercial  pesticides. 

Regardless  of  the  particular  construct,  the  average  proficiency  towards  paraoxon  is  3- 10-fold 
greater  than  towards  VX  (the  latter  being  based  on  the  DTNB  protocol).  The  best  variant  for 
paraoxon  hydrolysis,  114-2-Cl,  degraded  paraoxon  approximately  45-fold  faster  than  VX  (for 
VX  kinetic  data  see  our  6*  Year  Annual  Report).  This  variant  possesses  the  I117M  mutation, 
viz.,  of  a  residue  adjacent  to  the  active-site  residue,  W115,  and  in  a  proximity  that  may  influence 
the  binding  and  hydrolysis  of  paraoxon.  From  the  docking  model  of  paraoxon  in  the  active  site  of 
PONl,  it  seems  that  residue  PI  14  may  affect  the  positions  of  two  catalytically  important 
residues,  69  and  115,  and  that  residue  117  may  influence  the  positions  of  residues  114  and  115 
(Fig.  1). 


Figure  1.  Representation  of  paraoxon  docked  within  the  active  site  of  PONl  that  illustrates 
how  W115  and  adjacent  residues  may  influence  binding  of  paraoxon  and,  consequently,  its 
hydrolysis. 
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These  observations  call  for  the  reexamination  of  residues  114  and  117  in  the  search  for  an 
improved  catalyst  of  paraoxon  hydrolysis.  Specifically,  this  will  involve  the  generation  of 
libraries  that  will  explore  the  effects  of  substitutions  at  positions  1 14  and  1 17  on  the  efficiency  of 
paraoxon  hydrolysis. 

To  obtain  a  lead  variant  capable  of  effective  hydrolysis  of  malaoxon  and  phorate,  which  are 
defined  as  threats  to  the  civilian  population,  we  tested  the  capacity  to  hydrolyze  the  two 
pesticides  of  the  best  variants  in  both  the  rePON  1  and  rePTE  libraries  generated  for  hydrolysis  of 
G-  and  V-type  nerve  agents. 


S 

^2^50^  1 1 
C2H5O 


Phorate 


Although  malaoxon  contains  an  asymmetric  carbon  atom,  the  anti-AChE  activity  of  the  two 
enantiomers  differs  only  by  5-7-fold.  The  data  points  for  the  first  30-50%  of  hydrolysis  of 
malaoxon  by  rePONl  and  PTE  variants  fitted  mono-exponential  kinetics  very  well,  and  no 
attempts  were  made  to  resolve  the  individual  rates  for  each  enantiomer.  Activity  on  both 
pesticides  was  monitored  using  the  DTNB  protocol.  The  exact  concentration  of  each  pesticide 
was  determined  by  monitoring  the  total  release  of  the  corresponding  thiol  leaving  group  in  the 
presence  of  0.5  M  NaE/50  mM  phosphate,  pH  8.0.  Results  are  summarized  in  Tables  2  and  3  for 
rePON  and  rePTE  variants,  respectively. 

It  should  be  noted  that,  except  for  wt  PTE  and  PTE-C23,  due  to  the  extremely  slow  rates  of 
hydrolysis  the  reactions  were  recorded  up  to  30-50%  release  of  the  expected  amount  of  the  thiol 
leaving  group;  thus,  if  the  two  enantiomers  differ  in  their  rates  of  degradation  by  the  variants 
tested,  the  kj,^(/K„  values  reported  for  the  chiral  malaoxon  represent  the  faster  isomer. 
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Table  2:  Catalytic  activity  on  malaoxon  and  phorate  of  rePONl  variants  from  libraries 
VXG6  and  VXG7  (  xlO^  M'mm';  single  measurements). 


Library- Variant 

malaoxon 

phorate 

Wt  G3C9 

1.0 

<1  (no  reaction  at  2.1 
uM  over  60  min) 

VXG6  10H5 

5.1 

<2 

7F4 

2.0 

<2 

bG12 

2.5 

<2 

7F7 

4.5 

<2 

7C1 

1.7 

<2 

VXG7  167-F5 

8.3 

<2 

167-A6 

2.4 

<2 

223-C3 

2.0 

<2 

114-Cl 

6.7 

<2 

268-E8 

4.6 

<2 

Table  3:  Catalytic  activity  on  malaoxon  and  phorate  of  PTE 
(  xlO^  M  'min single  measurements). 


Library- Variant 

malaoxon 

phorate 

WtS5 

29  and  14 

(two  different  preps) 

19.5 

Early  libraries  119 

1.1 

C74 

2.1 

2.5 

A57 

8.1 

E36 

3.4 

3.3 

Library  G5  C23 

14 

14;  21 

A53 

8.2 

6.2 

Library  G6  IID9 

1.8 

2.0 

VF9 

4.9 

4.8 

IF4 

1.6 

4.6 

HID  10 

3.2 

4.4 

IVH3 

2.0 

5.5 

VA4 

2.6 

5.6 

variants 
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All  the  rePONs  tested  displayed  very  poor  activity,  if  any,  towards  phorate.  Under  the 
experimental  conditions  employed,  with  an  enzyme  protein  concentration  of  0.5-2  mM  in  the 
assay  cuvette,  the  activity  measured  was  indistinguishable  from  the  rate  of  spontaneous 
hydrolysis  of  the  combined  DTNB/phorate  blank  cuvette  (0.0024  min  ').  Thus,  we  estimated 
kcJKu  for  all  the  rePONs  listed  in  Table  2  to  be  <2x10^  M'mm'.  Phorate  is  probably  a  poor 
substrate  for  PONl  due  to  its  especially  unfavorable  structure,  which  contains  both  P=S  and  P-S- 
alkyl  moieties  (combining  both  thiono  and  thiolo  barriers).  In  the  case  of  malaoxon,  the  rate  of 
spontaneous  hydrolysis  of  DTNB/malaoxon  was  0.0028  min  ',  viz.,  only  2-3-fold  slower  than 
hydrolysis  in  the  presence  of  ~1  pM  rePONl,  thus  highlighting  the  fact  that  the  advanced 
rePONls  variants  tested  have  evolved  specifically  to  accommodate  the  CH3P(0)0-alkyl-  moiety 
of  VX,  and  are  less  suitable  for  binding  diaikyl-based  OPs  such  as  malaoxon  and  phorate. 

As  expected  from  the  catalytic  advantage  of  bacterial  PTEs  over  rePONs  when  reacting  with 
paraoxon  or  parathion,  the  PTE  variants  were  more  active  than  the  rePONs  towards  both  phorate 
and  malaoxon.  However,  the  observed  values  for  kj,^(/K„  listed  in  Table  3  indicate  that  the 
variations  introduced  into  the  PTE  scaffold  so  as  to  accelerate  the  hydrolysis  of  V-agents,  did  not 
allow  for  sufficient  catalytic  activity  on  either  phorate  or  malaoxon  to  permit  their  consideration 
as  antidote  candidates  for  either  of  these  OPS.  In  fact,  with  the  exception  of  C23,  in  all  cases  the 
mutations  introduced  decreased  catalytic  activity  relative  to  the  wt  enzyme.  It  should  be  noted 
that,  despite  the  potent  catalysis  of  paraoxon  and  parathion  (both  diethyl  phosphate-based 
pesticides)  by  the  wt  PTE  variant  S5,  its  activity  towards  the  dime  thy  Iphosphoryl  analog  is 
reduced  significantly,  presumably  due  to  lack  of  tight  binding  relative  to  the  diethy Iphosphoryl 
homolog.  Thus,  we  found  that  parathion  (containing  diethoxyphosphoryl)  is  hydrolyzed  by  wt  S5 
20-fold  faster  than  its  dimethoxy  homolog,  methyl  parathion;  this  despite  the  greater  electron 
donating  properties  of  the  ethoxy  relative  to  the  methoxy  group,  which  should  decrease  the 
electropositivity  of  the  P  atom,  thus  making  it  an  inferior  substrate  to  methyl  parathion.  Indeed, 
fluoride  was  found  to  hydrolyze  parathion  3-fold  more  slowly  than  methyl  parathion.  Thus,  the 
low  kj,^(/K„  values  observed  for  hydrolysis  of  malaoxon  by  the  PTEs  are  attributed  to  the 
combined  effects  of  poor  molecular  complementarity  of  the  malaoxon  molecule  and  of  the  thiolo 
barrier. 
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Methyl  parathion  parathion 

As  expected,  the  conversion  of  the  P=0  of  malaoxon  to  the  thiono  function  (P=S)  of 
malathion  further  decreased  the  susceptibility  of  this  OP  type  to  both  rePTEs  and  rePONs. 


malathion  malaoxon 


Thus,  when  PTE  variants  of  library  6  (IID9,  VE9,  IE4,  HID  10,  IVH3,  and  VA4,  at  0.4- 1.0 
pM)  were  reacted  with  malathion,  the  activity  measured  was  indistinguishable  from  the 
spontaneous  rate  of  chromophore  release  in  the  combined  DTNB /malathion  control  run  (0.0015 
min  ');  we  therefore  estimate  that  of  these  variants  is  <2x10^  M  'min  '.  The  five  best 

rePONl  variants  of  library  VXG7  (114-Cl,  167-A6,  167-E5,  223-C3,  and  268-E8)  also  did  not 
display  a  measurable  increased  activity  over  control  values  towards  malathion,  despite  the  fact 
that  they  were  at  1.5-3  pM  in  the  final  reaction  mixture.  We  conclude  that  the  dimethoxy 
phosphoryl  pesticides  are  very  poor  substrates  of  the  OP  hydrolases  examined,  and  that  a 
different  line  of  study  is  required  to  evolve  variants  with  sufficient  activity  to  qualify  as 
bioscavengers  of  the  dimethylphosphoryl-based  OPs. 

#2.  Validation  of  the  hypothesis  concerning  the  structural  motifs  in  PONl  that 
are  involved  in  association  with  HDL,  improving  rePONl ’s  affinity  for  HD, 
and  increasing  the  biological  life  time  of  rePON  variants 

In  order  to  increase  the  stability  of  rePONl  variants  in  the  circulation  we  examined  the 
influence  of  the  known  capacity  of  PONl  to  interact  with  HDE,  which  is  accompanied  by  an 
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increase  in  its  catalytic  activity.  We  earlier  hypothesized  that  the  H2  helix  and  its  connecting 
loops,  the  HI  helix,  and  the  flexible  loop  residues,  are  together  responsible  for  the  structural 
changes  that  result  in  stimulation  of  activity  when  PON  1  is  bound  to  HDL  (see  our  2"‘^  Annual 


Report). 

In  order  to  test  the  importance  of  individual  residues  on  HDL  binding  we  developed  a  cell- 
free  system  that  permits  the  synthesis  of  rePON  1  and  its  various  mutants  in  the  presence  of  HDL 
(depicted  in  the  cartoon  below).  The  advantage  of  this  technique  is  that  it  alleviates  both  the  need 
for  protein  purification  in  detergent  buffers,  and  for  subsequent  removal  of  detergent  prior  to 
testing  the  interaction  of  variants  with  HDL.  Indeed,  removal  of  detergent  from  purified  rePONl 
was  found  to  greatly  compromise  its  stability.  In  addition,  this  assay  system  can  be  used  to 
analyze  many  mutants  in  parallel.  The  selection  of  putative  key  positions  in  PONl’s  helices  was 


based  on  structural  analysis  of  our  rePONl  crystal  structures,  as  well  as  on  sequence  analysis  of 


the  PON  family. 


Reaction  mix 

„  .  Enzymes,  amino-acids, 

Reconstituted  ,  „  ,  , 

high-density  iipoprotein  (HDL)  ^  ^ 


V 

Ceii-free  expression  in  batch  mode: 

•  Volume:  15-30^1 

•  Expression-time:  2h 

•  Expression  yield:  “100-500pg/ml 


E.  Coli  Extract  (S30) 


Recombinant  piasmid 


In  our  6*  year  Annual  Report  we  summarized  the  preliminary  data  that  pointed  to  the 
involvement  of  Y185  and  F186,  within  the  H2  connecting  loop,  especially  of  F186,  which 
displays  a  significant  decrease  in  the  stimulation  of  activity  towards  TBBL.  We  also  pointed  out 
that  residue  Y71,  within  the  flexible  loop,  which  moves  significantly  upon  binding  of  2HQ,  is 
likely  to  be  a  key  residue  in  stimulation  of  the  lipo-lactonase  activity  of  rePONl.  Furthermore, 
unlike  the  L198Y  mutation,  which  results  in  both  stabilization  and  stimulation  of  activity  by  this 
H2  residue,  the  L198W  mutation  results  only  in  stimulation  of  activity  relative  to  the  WT. 
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We  report  here  on  the  evaluation  of  positions  in  the  HI  helix  and  in  the  active-site  region  of 
the  wt  G3C9  variant,  for  which  results  are  summarized  in  Table  4. 


Table  4:  Mutations  in  helix  HI  and  within  the  active  site  produce  enhanced  stability  and 
activity  relative  to  WT  rePON  1  in  the  presence  of  HDL 


Variant 

Activity  stimulation 

Inactivation  144hr 

^app(nM) 

(%) 

Residual  lactonase 
activity  (%) 

rePONl 

0.58±0.07 

406±11 

61±6 

A6V 

0.53±0.07 

505±15 

84±3 

LTV 

1.3±0.17 

364±13 

42±5 

T8V 

0.41±0.05 

508±11 

73±7 

T8I 

0.44±0.01 

450±3 

66±4 

G13V 

0.58±0.05 

500±12 

74±5 

Y71F 

0.38±0.1 

467±14 

41±4 

Y71G 

1.70±0.09 

430±7 

35±6 

H184Q 

1.5±0.2 

367±12 

6±2 

H184T 

1.00±0.09 

613±14 

84±3 

The  first  row  entries,  for  rePONl,  are  relative  to  the  same  amount  of  protein  in  the  absence  of  HDL. 
Activity  stimulation  refers  to  TBBL  hydrolysis,  and  inactivation  tests  were  performed  in  the  presence  of 
15  mM  EDTA  at  37°C.  Residues  6,  7,  8  and  13  are  on  the  HI  helix,  and  71  and  184  are  located  within  the 
active  site. 

As  can  be  seen,  stimulation  of  activity  for  wt  rePONl  in  the  presence  of  HDL  was  ~4-fold. 
Mutations  at  positions  6  and  13  of  Hi,  and  at  active-site  residue  184,  further  increased  activity. 
This  enhancement  was  accompanied  by  increased  stability  for  mutations  at  positions  6,  13  and 
184. 

Thus  far,  our  conclusions  were  based  on  the  stabilization  and  enhanced  activity  revealed  by 
monitoring  lactonase  activity.  To  demonstrate  the  broad-spectrum  enhancement  achieved  by 
association  with  HDL,  the  activity  of  several  variants  obtained  by  mutating  the  G3C9  scaffold 
was  monitored  towards  the  VX  surrogate,  BMP,  in  which  the  diisopropylamino-ethanthiolo 
leaving  group  had  been  replaced  by  coumarin.  As  shown  in  Fig.  2,  increasing  the  HDL 
concentration  produced  complexes  of  both  the  G13V  and  H184T  mutants  that  displayed 
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increased  activity  relative  to  the  wt  enzyme.  Deletion  of  the  20  amino  acids  at  the  N-terminus 
decreased  activity,  and  thus  substantiated  the  importance  of  this  region  in  binding  to  HDL 


□ 


0.5n 


HDL,  pM 


wtG3C9 

L198A 

delta  20  (N-terminus) 

G13V 

H184T 


Figure  2:  Increased  binding  of  PONl  to  HDL  increases  its  catalytic  activity  on  the  VX 
surrogate,  BMP. 


When  the  best  mutated  positions  in  wt  G3C9  were  copied  onto  the  IIGl  scaffold  (which  is 
the  most  effective  variant  towards  cyclosarin  and  soman),  and  tested  on  the  cyclosarin  surrogate, 
CMP,  we  obtained  essentially  the  same  pattern  of  enhancement,  albeit  with  a  reduced  efficacy 
compared  to  the  wt  G3C9  (see  Fig.  3  below).  Thus,  the  effects  of  modification  of  the  HI  and  H2 
helices,  together  with  the  connecting  loops,  are  dependent  on  other  residues  yet  to  be  identified. 
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Figure  3:  Stimulation  of  hydrolysis  of  the  GF  analog  CMP  due  to  binding  to  HDL 


In  order  to  examine  whether  mutations  that  can  stimulate  rePONl  activity  and  enhance  its 
stability  in-vitro  can  also  provide  it  with  an  increased  circulatory  half-life  in-vivo,  we  sent  the 
following  variants  for  large-scale  production  and  purification  by  the  Magliery  lab  at  OSU,  to  be 
followed  by  in-vivo  pharmacokinetics  testing  by  the  Cerasoli  lab  at  ICD:  A6V-G3C9  and  G13V- 
G3C9,  variants  that  display  increased  HDL  binding,  and  Del20-G3C9,  a  variant  that  displays 
decreased  HDL  binding  relative  to  G3C9,  and  thus  provides  a  negative  control.  By  comparing 
the  in-vivo  circulatory  half-life  times  of  these  proteins  relative  to  wt-rePONl  (G3C9),  we  should 
be  able  to  validate  our  hypothesis,  and  provide  an  approach  to  enhancing  the  pharmacokinetic 
properties  of  rePONl  variants. 

On  August  26*,  2013,  all  three  rePONl  clones  (A6V-G3C9,  G13V-G3C9  and  N-Del20- 
G3C9)  were  sent  to  OSU  for  expression  and  purification,  with  the  aim  of  studying  their  PK 
profile  in  guinea  pigs  at  ICD.  We  are  awaiting  the  results  of  the  in  vivo  experiments  so  as  to 
design  the  next  generation  of  variants  with  increased  affinity  for  HDL. 
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In  the  appendix  we  summarize  a  comprehensive  study  that  describes  the  structural  and 
mechanistic  basis  for  the  HDL- induced  stimulation  of  PONl’s  activity. 

#3.  Improved  modeling  of  VX  within  the  rePONl  scaffold 

The  objective  of  this  modification  was  to  enhance  the  accuracy  of  Steered  Molecular 
Dynamics  (SMD)  and  Molecular  Dynamics  (MD)  simulations  for  docking  of  VX,  in  order  to 
achieve  improved  prediction  of  the  properties  of  newly  designed  rePONl  mutants.  Fig.  14 
models  the  complex  of  1-6-G4  (from  round  VXG4)  with  the  toxic  Sp-VX  enantiomer.  This 
variant  is  ~  100-fold  more  active  toward  VX  than  the  wt  enzyme.  The  mutations  in  the  variant 
were  copied  into  the  G2E6  crystal  structure,  with  the  exception  of  two  residues  that  are  not 
within  the  active  site: 

□ 

1-6-G4;  L55I,  L69V,  K70N,  Y71M,  H115W,  H134R,  H184A,  K192Q,  F222M,T332S 


The  model  was  generated  with  a  new  simulation  approach,  and  the  most  important 
differences  from  our  previous  modeling  are  as  follows: 

1)  The  protein  and  ligand  were  parameterized  by  using  Amber99sb-ildn  force  field  [1]  rather 
than  the  OPLS/AA  force  field  [2,3]  that  was  used  previously.  It  has  been  shown  that  simulations 
of  proteins  parameterized  with  the  Amber99sb-ildn  force  field  reproduce  results  of  NMR  studies 
better  than  any  other  modern  force  field  [4] . 

2)  Another  advantage  of  the  Amber99sb-ildn  force  field  is  the  possibility  of  introducing  a 
more  sophisticated  model  of  water  -  TIP5P  [5],  which  is  not  possible  with  OPLS/AA.  TIP5P 
reproduces  the  electrostatic  properties  of  water  better  than  the  previously  used  TIP4P  [6].  In 
TIP5P,  lone  electron  pairs  on  the  water  oxygen  are  represented  by  pseudo-atoms.  This  permits 
improved  accuracy  of  simulations  of  biomolecular  interaction  phenomena  in  which  microscopic 
behavior  of  water  molecules  plays  an  important  role,  as  is  the  case  in  PONl-mediated  catalytic 
OP  degradation  (Fig.  1). 

3)  A  new  simulating  pipeline  was  deployed  on  a  fast  GPU  device,  which  increased  the  speed 
of  simulation  significantly. 
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It  is  envisaged  that  the  overall  increase  in  the  reliability  of  the  simulation  methodology  will 
improve  the  accuracy  of  the  prediction  of  the  catalytic  properties  of  new,  rationally  designed 
variants  of  rePONls  when  reacting  with  V-type  nerve  agents.  Compared  to  previous  attempts, 
the  current  results  generate  a  more  detailed  description  of  the  initiation  of  the  catalytic  process 
and  provide  better  input  for  further  study,  for  example  with  QM/MM  methods. 


Figure  14:  Model  of  the  1-6-G4  rePONl  variant  complexed  with  Sp-VX. 


The  orientation  of  VX  fulfills  the  criteria  of  the  binding  mode  necessary  for  catalytic 
breakdown  of  OP  according  to  the  mechanism  that  is  considered  the  most  probable  [7] .  In  this 
mechanism,  a  hydroxyl  ion  produced  by  deprotonation  of  water  "1"  by  D269  and  E53  (lines  "a" 
and  "b")  attacks  the  phosphorus  atom  of  VX  (electrostatically  attracted  to  the  catalytic  calcium 
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ion,  Ca357  -  line  "c")  along  the  line  "d".  As  a  consequence,  the  S-P  bond,  marked  by  the  arrow, 
is  severed.  The  leaving  group  is  stabilized  by  a  possible  H-bond  with  N168.  Structural  water  "2" 
is  conserved  in  all  the  rePONl  crystal  structures.  It  is  proposed  that  this  water  molecule  mediates 
transfer  of  the  proton  of  the  carboxyl  group  of  residue  E53,  and  restores  its  structure  in  the  native 
enzyme  so  as  to  initiate  a  new  catalytic  cycle. 

Key  Research  Accomplishments 

•  Identification  of  a  variant  catalyzing  hydrolysis  of  paraoxon  at  2x10^  M  'min  ' 

•  Elucidation  of  the  structural  and  mechanistic  basis  for  the  HDE-induced  stimulation  of 
PONl’s  lactonase  activity.  Improved  capacity  to  associate  rePONIs  with  HDE, 

•  Use  of  steered  molecular  dynamics  (SMD)  methodology  to  model  the  interaction  of  VX 
with  PONl 

Conclusions 

•  Screening  existing  libraries  with  paraoxon,  parathion,  malaoxon,  malathion,  and  phorate 
showed  that  with  the  exception  of  paraoxon,  all  the  rePONs  tested  displayed  poor 
activity,  if  any,  and  a  different  line  of  enhanced  evolution  will  have  to  performed  so  as  to 
generate  variants  capable  of  hydrolyzing  them  effectively. 

•  The  improved  modeling  of  VX  into  rePON  1 ,  so  as  to  achieve  enhanced  proficiency  in  the 
design  of  new  variants  suggests  several  sites  for  introducing  mutations  that  will  enhance 
rates  of  degradation  of  V-type  nerve  agents. 

•  The  comprehensive  study  of  rePONl -HDE  association  identified  residues  involved  in 
both  enhancement  of  catalysis,  as  well  as  residues  that  produced  enhanced  thermal 
stability,  offering  new  approaches  to  extension  of  the  biological  life-time  of  rePONs. 
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Appendix-Supplementary 


Model  for  the  binding  of  PONl  to  HDL 

In  this  appendix  we  provide  a  detailed  description  of  a  model  for  the  binding  of  PONl  to 
HDL,  especially  with  regard  to  the  H2  residues,  obtained  via  mutational  and  structural  analyses 
of  PONl  along  with  computational  simulations.  We  thus  describe  the  structural  and  mechanistic 
basis  for  the  HDL-induced  stimulation  of  PONl’s  lactonase  activity.  Furthermore,  we  discuss 
this  model  within  the  framework  of  a  more  general  hypothesis  regarding  the  importance  of 
structural  rigidity,  and  the  importance  of  second-  and  third-shell  residues,  in  configuring  enzyme 
active  sites. 

The  PONl/HDL  structural  model 

We  already  earlier  postulated  the  interactions  of  PONl’s  helices  HI  and  H2  with  the  surface 
of  HDL,  based  primarily  on  the  presence  of  exposed,  surface  hydrophobic  residues,  especially  of 
aromatic  amino  acids  (Gaidukov  &  Tawfik,  2005;  Harel  et  al,  2004)  (Figure  SI).  However,  the 
proposed  model  did  not  address  the  key  question  -  namely,  how  does  HDL  binding  induce 
dramatically  higher  stability  and  catalytic  efficiency.  Our  working  hypothesis  has  been  that  HDL 
binding  would  immobilize  the  PONl  segments  attached  to  it,  which,  in  turn,  would  affect  the 
active-site’s  shape  and/or  catalytic  organization.  In  this  context,  we  have  now  examined  the  B- 
factors  of  available  PONl  structures,  and  found  that  the  C-terminal  end  of  HI  (the  N-terminus  is 
completely  disordered)  and  H2  comprise  the  most  mobile  parts  of  the  enzyme.  In  accordance,  the 
normal  mode  analysis  of  these  structures  (NMA;  ElNemo  server  -  http://igs-server.cnrs- 
mrs.fr/elnemo/)  revealed  that  helices  HI  and  H2  show  large,  coordinated  motions.  This  motion 
results  in  very  different  configurations  of  the  active  site’s  opening,  and  in  the  displacement  of 
catalytic  residues,  most  notably  N168  (a  Ca^^- ligating  residue)  and  the  catalytic  Ca^^  (Figure 
S2A-B). 

Further  analysis  of  the  crystal  structures  of  rePONls  in  the  region  of  helices  HI  and  H2 
revealed  potential  interactions  of  a  molecule  of  the  detergent  n-dodecyl-beta-D-maltoside 
(DDM)  that  was  present  in  the  crystallization  medium.  The  detergent  seems  to  interact  with 
residues  of  HI  (e.g.,  Q20),  as  well  as  with  the  adjacent  H2’s  N-terminal  connecting  loop  (L200 
and  W202;  Figure  S2C).  The  detergent,  and  presumably  the  lipid  components  of  HDL,  may  thus 
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restrict  these  H1-H2  motions.  Indeed,  the  packing  surfaces  between  rePONl  molecules  within 
the  different  asymmetric  units  of  the  unit  cell  involve  close  interactions  between  the  adjacent  H2 
helices,  as  well  as  between  their  connecting  loops,  and  specifically  between  residues  Y185  and 
F186  of  the  neighboring  PONl  molecules  (Figure  S2D).  It  appears,  therefore,  that  both  the 
bound  DDM  molecules,  and  the  crystal  packing  interactions  involving  HI  and  H2,  play  a  major 
role  in  limiting  PONl’s  conformational  flexibility  and,  appear  to  affect  the  way  crystals  form. 
Given  the  hydrophobic  nature  of  these  interactions,  they  are  likely  to  mimic  the  binding 
interactions  of  PONl  with  HDL.  Furthermore,  the  HI  and  H2  configurations  seen  in  the  crystal 
structure  may  represent  only  one  out  of  many  possible  configurations  available  for  these  helices, 
and  hence  of  PONl’s  active-site.  Since  conformational  flexibility  promotes  instability  and  faster 
unfolding  rates  (Matthews  et  al,  1987),  HDL  binding  increases  the  enzyme’s  thermodynamic 
stability  (as  does  detergent  binding,  though  to  a  lesser  degree).  Similarly,  the  preferred  adoption 
of  a  specific  active-site  configuration  may  account  for  the  HDL-induced  increase  in  catalytic 
efficiency. 

In  situ  incorporation  of  rePONl  into  reconstituted  HDL  particles 

The  experimental  protocol  applied  thus  far  for  the  analysis  of  PONl’s  interaction  with  HDL 
includes:  1)  production  and/or  purification  of  recombinant  or  of  human  PONl  in  the  presence  of 
a  detergent;  2)  delipidation  for  removal  of  the  detergent;  3)  incubation  of  the  detergent-depleted 
PONl  with  HDL.  Aside  from  being  laborious,  this  protocol  is  problematic  due  to  the  rapid 
inactivation  of  the  delipidated  PONl.  Thus,  the  observed  catalytic  stimulation  might  be  over¬ 
estimated  when  compared  to  a  partly  inactive  control  devoid  of  lipids.  To  enable  faster  and  more 
reliable  analysis,  we  developed  a  cell-free  expression  system  that  permits  the  expression  of 
rePONl  in  the  presence  of  reconstituted  HDL  particles  (rHDL).  The  latter  were  prepared  as 
discoidal  rHDL  particles  containing  l-palmitoyl-2-oleoylphosphatidylcholine  (POPC),  free 
cholesterol  (FC),  and  the  apolipoprotein  apoA-I,  at  a  molar  ratio  of  100/5/1  (Gaidukov  et  al, 
2010;  Oda  et  al,  2001).  We  first  validated  the  reliability  of  this  procedure  by  testing  the 
stimulation  of  activity  of  the  wt-like  rePONl  towards  a  chromogenic  lactone  substrate,  TBBL, 
and  towards  paraoxon.  We  measured  the  activity  of  the  cell-free  expression  products  at  0-10  pM 
rHDL. 
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The  results  obtained  with  this  novel  procedure  were  in  agreement  with  those  obtained  with 
the  previously  applied  delipidation  protocol  (Gaidukov  &  Tawfik,  2005;  Gaidukov  et  al,  2010). 
In  essence,  PON  1 -catalyzed  TBBL  hydrolysis  was  stimulated  by  HDL,  following  a  typical 
Langmuir  isotherm,  and  the  values  (percentage  activity  relative  to  the  activity  of  the  cell- 
free  translated  PON  1  in  the  absence  of  rHDL)  and  (the  rHDL  concentration  that  produces 
50%  of  the  total  stimulation  of  activity)  displayed  similar  trends.  As  previously  observed,  the 
paraoxonase  activity  showed  only  minor  stimulation,  and  stimulation  of  lactonase  activity  was 
completely  abolished  in  the  HI -truncated  form. 

The  effect  of  rHDL  on  the  stability  of  rePONl  was  explored  by  a  previously  applied 
inactivation  assay.  The  cell-free  expressed  PONl  samples  in  the  absence  of  HDL,  and  at  the 
highest  rHDL  concentration  (lOpM),  were  incubated  at  37°C  in  the  presence  of  a  calcium 
chelator,  EDTA.  The  residual  lactonase  activity  was  tested  after  incubation  with  EDTA  for  72  h, 
and  compared  to  the  activity  prior  to  EDTA  addition  (Table  SI).  The  decreased  inactivation 
rates  upon  HDE  binding  (Gaidukov  &  Tawfik,  2005;  Gaidukov  et  al,  2010)  were  similarly 
reproduced.  Table  SI  reveals  that  lactonase  stimulation  and  stabilization  are  overall  correlated. 

Mutational  analysis  of  H2  and  of  its  C -terminal  loop 

Truncation  of  Hi  resulted  in  a  ~100-fold  drop  in  affinity  for  HDE,  and  in  abolition  of  both 
HDE- mediated  stabilization  and  stimulation  of  enzymic  activity,  thus  providing  evidence  for  the 
role  of  Hi  in  the  HDE  binding  (Gaidukov  &  Tawfik,  2005)  (Table  SI).  Since  complete 
truncation  of  H2  is  not  an  option,  point  mutations  were  explored  with  the  aim  of  examining  the 
role  of  H2  in  HDE  binding. 

H2  and  its  C-terminus  loop  possess  a  stretch  of  hydrophobic  surface  side-chains  that  were 
postulated  to  be  involved  in  the  binding  to  HDE  (Gaidukov  &  Tawfik,  2005;  Harel  et  al,  2004). 
These  include  H2  residues  Y190,  W194  and  E198,  and  the  C-terminal  loop  residues  Y185  and 
El 86  (Figure  SI).  The  mammalian  PONs  exhibit  high  sequence  conservation  at  these  positions  - 
a  degree  of  conservation  that  is  unexpected  for  surface  residues.  Some  positions  are  fully 
conserved  {e.g.,  positions  183-186),  and  others  display  only  conservative  substitutions  {e.g., 
position  190,  Tyr  vs.  Phe).  Polio  wing  the  structural  and  sequence  analyses,  we  designed  mutants 
that  sample  H2  residues  Y190,  W194  and  E198,  as  well  as  Y185  and  P186  within  H2’s  C- 
terminal  loop  (Figure  SI).  Substitutions  to  Ala,  especially  of  positions  190  and  194,  reduced  the 
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HDL- induced  stabilization  and  lactonase  stimulation  (Table  1),  with  the  190/194  double  mutant 
showing  almost  no  stimulation  or  stabilization  {V„aJKapp  values  reduced  by  ~14-fold).  Mutations 
in  H2’s  N-terminal  loop  also  had  a  marked  effect.  While  Y185A  showed  comparable  parameters 
to  that  of  the  wild  type,  Y185G,  that  is  likely  to  induce  changes  not  only  in  the  side-chain,  but 
also  in  the  loop’s  backbone  configuration,  resulted  in  near-complete  loss  of  the  effects  of  HDL. 
Both  F186A  and  F186G  resulted  in  a  decrease  in  lactonase  stimulation  and  stability  (Table  SI). 
The  double  mutant,  Y185A/F186A,  showed  a  >10-fold  decrease  in  both  lactonase  stimulation 
and  stability  (Table  SI).  Overall,  these  results  support  the  hypothesis  that  residues  of  helix  H2, 
as  well  as  of  its  connecting  loop,  play  a  role  in  HDL  binding. 

The  mechanism  of  PONl’s  catalytic  stimulation 

PONl’s  crystal  structures  reveal  a  long  loop  that  connects  the  N-terminus  of  H2  with  the 
enzyme’s  scaffold  and  active  site.  This  loop  also  includes  the  active-site  residues  H184  and 
D183  (Figures  S2B  and  S3).  Mutations  in  these  two  residues  render  both  human  PONl  and 
rePONl  almost  inactive  (Josse  et  al,  1999;  Khersonsky  &  Tawfik,  2006).  However,  a  detailed 
analysis  indicated  that  neither  of  these  residues  is  directly  involved  in  catalysis,  i.e.,  as  a  base  or 
acid  catalyst  (Khersonsky  &  Tawfik,  2006).  Indeed,  analysis  of  the  crystal  structure  of  the  wild- 
type-like  rePONl  in  complex  with  a  lactone  analogue  revealed  no  direct  contacts  of  the  substrate 
with  either  D183  or  H184,  and  hence  the  proposed  catalytic  mechanism  did  not  include  these 
residues  (Ben-David  et  al,  2012).  These  residues,  especially  D183,  are  highly  conserved,  and  are 
found  in  very  distant  members  of  the  PON  family  {e.g.,  bacterial  PONs  with  <30%  sequence 
similarity  to  mammalian  PONs)  (Bar-Rogovsky  et  al,  2013).  What  then  are  the  roles  of  H184 
and D 183? 

Despite  the  absence  of  direct  contacts,  there  exists  a  network  of  hydrogen  bonds  spanning 
from  H2’s  C-terminus  loop,  and  particularly  residues  H184  and  D183,  to  two  active-site  residues 
that  are  in  direct  contact  with  the  substrate:  residue  Y71  -  at  the  tip  of  a  highly  mobile  active-site 
loop  (residues  70-81),  and  the  catalytic  Ca^^-ligating  residue,  N168  (Figure  3S).  As  shown  (Ben- 
David  et  al,  2012;  Ben-David  et  al,  2013),  conformational  variations  in  the  mobile  active-site 
loop,  as  well  as  in  the  catalytic  Ca^^,  modulate  rePONl ’s  catalytic  activity,  and  specifically  relate 
to  the  enzyme’s  dual  activity  as  a  lactonase  (the  native  function)  and  as  a  paraoxonase  (a 
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promiscuous  activity).  We  thus  reasoned  that  the  lactonase  stimulation  seen  upon  HDL  binding 
might  be  mediated  via  the  hydrogen  bonds  of  D 183  to  Y71  and/or  to  N 168  (Figure  S3). 

To  examine  these  two  options,  we  prepared  a  set  of  mutants  in  residues  that  are  part  of  the 
network  that  includes  D183-H184,  Y71  and  N168  (Figure  S3).  As  previously  reported,  all 
mutations  at  positions  183  and  168  resulted  in  complete  loss  of  both  the  lactonase  and 
paraoxonase  activities.  Other  mutants  did  show  detectable  activity,  including  Y71F,  Y71G, 
H184Q  and  H184T.  These  mutations  had  a  relatively  mild  effect  on  the  lactonase  activity  (<6- 
fold  decrease  towards  TBBL  (Khersonsky  &  Tawfik,  2006)).  We  thus  tested  the  HDL- induced 
stimulation  of  the  activity  of  these  mutants.  The  conservative  Y71F,  but  also  the  drastic  Y71G 
mutation,  resulted  in  only  a  mild  decrease  in  the  HDL-induced  stabilization  and  catalytic 
stimulation  (<3-fold;  Table  SI).  Foremost,  abolishing  the  putative  H-bond  between  Y71’s 
hydroxyl  (via  the  Y71F  mutation)  and  D183  (Figure  3)  had  no  effect  on  the  HDL-induced 
stimulation  of  activity  (Table  SI).  This  suggests  that  the  putative  hydrogen  bonds  of  183D-184H 
with  Y71  do  not  play  a  key  role  in  PONl’s  activity  in  general,  and  specifically  in  the  HDL- 
induced  stimulation. 

The  hydrogen  bonds  spanning  from  H2,  including  K192,  to  N168,  are  much  more  likely  to  be 
crucial  for  catalysis.  This  hypothesis  is  supported  by  the  absolute  essentiality  of  N 168  (mutants 
of  other  residues  that  ligate  the  Ca^^  -  53,  224,  269  and  270  -  do  display  some  residual  activity) 
(Ben-David  et  al,  2013))  and  of  D183  (Khersonsky  &  Tawfik,  2006),  and  by  the  known  effects 
of  the  192Q/R  polymorphism  on  HDL  binding  and  stimulation  of  activity. 

The  role  of  HI 84  in  aligning  D183,  and  in  turn  N 168,  is  supported  by  the  fact  that  mutations 
of  H184  to  non-polar  residues,  e.g.,  to  Ala  (Josse  et  al,  1999;  Khersonsky  &  Tawfik,  2006), 
result  in  loss  of  activity,  but  polar  residues  such  as  glutamine  or  threonine  retain  activity  and 
stimulation,  or  even  slightly  increase  the  stimulatory  effect  of  HDL  (H184T;  Table  SI).  Further, 
PONl’s  denaturation  is  initiated  by  the  loss  of  the  catalytic  Ca^^,  followed  by  loss  of  the 
structural  one,  and  irreversible  unfolding  (Gaidukov  &  Tawfik,  2005;  Kuo  &  La  Du,  1998).  The 
spatial  fixation  of  N 168  by  D183,  and  the  accompanying  electrostatic  effects  that  may  stabilize 
the  catalytic  Ca^^,  make  a  significant  contribution  to  the  affinity  of  the  catalytic  Ca^^for  PONl. 
Thus,  limiting  the  conformational  freedom  at  the  H1-H2  interface  (Figure  S2),  and  thereby 
reinforcing  the  H-bond  network  that  immobilizes  N168  and,  in  turn,  the  catalytic  Ca^^  (Figure 
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S3),  provides  a  likely  explanation  for  the  effects  of  HDL-binding  on  both  PONl’s  lactonase 
activity  and  its  configurational  stability. 

EVB  simulations 

A  crystal  structure  of  PONl  bound  to  an  HDL  particle,  or  for  that  matter,  of  any  other  HDL- 
associated  enzyme,  is  unavailable.  Thus,  the  structural  changes  induced  by  HDL  binding,  and 
their  effect  on  the  active-site  configuration,  remain  unknown.  To  gain  further  insights  regarding 
the  mode  of  stimulation  of  PONl  activity,  we  therefore  used  the  crystal  structure  of  the  wild- 
type-like  rePONl  as  the  basis  for  empirical  valence  bond  (EVB)  simulations.  Our  working 
assumption  has  been  that  the  changes  in  PONl’s  lactonase  activity  upon  HDL  binding  relate 
primarily  to  the  conformation  of  H2,  as  indicated  by  the  experimental  data  presented  above.  A 
series  of  EBV  models  were  generated  by:  (1)  simulating  the  gradual  pulling  of  H2  away  from  the 
active  site;  (2)  using  a  range  of  H2  configurations  obtained  from  normal  mode  simulations 
(Figure  SIA,  B).  Simulations  were  performed  in  the  presence  of  a  lactone  substrate  (TBBL)  that 
represents  PONl’s  lipophilic  lactone  substrates,  as  well  as  with  paraoxon,  a  promiscuous 
substrate.  The  models  were  simulated  and  relaxed  for  100  ps  at  300K. 

Two  alternative  models  for  the  TBBL/PONl  complex  were  examined.  Both  models  were 
based  on  HI  15  as  the  general  base  generating  the  hydroxide  nucleophile  (Ben-David  et  al,  2012). 
The  models  differed  by  either  N168  or  H285  providing  leaving  group  assistance  (LGA),  and  the 
substrate’s  orientation  was  refined  accordingly.  We  found  that  during  relaxation,  the  structures 
with  either  no  perturbation  or  minor  perturbation  of  H2’s  position  maintained  well-ordered  active 
sites.  Specifically,  N168  and  N224  maintained  a  close  interaction  with  the  both  the  metal  ion  and 
the  substrate.  However,  upon  more  significant  perturbation  of  H2,  the  side-chain  positions  of 
both  H184  and  D183  were  disturbed,  and  subsequently,  the  catalytic  Ca^^’s  interaction  with 
N168  was  severely  altered  (Figure  S4).  Simulations  were  also  performed  with  H2  constrained 
throughout  the  relaxation  step  and  EBV  runs.  This  might  suggest  that  H285  only  serves  as  an 
auxiliary,  backup  residue  (Ben-David  et  al,  2012)  in  cases  in  which  N168  is  highly  fluctuating 
(i.e.,  when  PONl  is  not  bound  to  HDE). 

EVB  models  with  paraoxon  were  also  examined,  with  D269  as  the  general  base  (Ben-David 
et  al,  2013).  As  seen  with  the  TBBE  models,  perturbing  the  position  of  H2  affected  the  network 
of  interactions  that  includes  H184,  D183  and  N168.  Nonetheless,  the  active  site  maintained  its 
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catalytic  configuration  for  paraoxon  hydrolysis.  Notably,  N224’s  side-chain  showed  some 
flexibility,  whereas,  as  opposed  to  the  lactone  hydrolysis  trajectories,  N168  largely  remained 
fixed.  The  fixation  of  N168  despite  H2’s  conformational  freedom  seems  to  be  driven  by  the 
interaction  of  N168  with  paraoxon’s  well-positioned  leaving  group  oxygen.  Thus,  paraoxon 
hydrolysis  seems  to  be  less  sensitive  to  the  active  site’s  configurational  ensemble,  which  is  in 
line  with  paraoxon  being  a  promiscuous  substrate,  and  with  HDL  binding  inducing  a  minor 
change  in  PONl’s  paraoxonase  activity  -  ~1. 5-fold  stimulation,  as  compared  to  <20-fold  towards 
lactones  (Gaidukov  &  Tawfik,  2005). 

Crystal  structure  of  the  192Q  mutant 

The  192RyQ  variants  are  polymorphs  of  human  PONl  that  were  reported  to  differ  in  their 
antiatherogenic  potential,  with  the  192R  isozyme  displaying  higher  potency  (Costa  et  al,  2013; 
She  et  al,  2012).  It  was  accordingly  found  that  the  192Q  polymorph  exhibits  significantly 
reduced  HDL  binding  and  lactonase  stimulation.  Similarly  to  rabbit  PONl,  rePONl  has  Lys  at 
position  192.  RePONl-192K  and  -192R  show  essentially  the  same  level  of  stability  and 
activation  upon  the  association  to  HDL,  whereas,  similarly  to  the  human  polymorph,  re-PONl- 
192Q  shows  reduced  stability  and  stimulation  (Gaidukov  et  al,  2006).  However,  the  structural 
basis  for  the  different  behavior  of  192K/R  versus  192Q  remained  thus  far  unexplained. 

Previously  solved  crystal  structures  of  rePONl-192K  indicate  that  K192  resides  on  H2,  with 
the  side-chain  pointing  towards  the  active  site  (Figure  S3;  PDB  codes:  XXXX).  We  have  now 
determined  the  crystal  structure  of  re-PON  1-192Q.  As  in  our  previous  studies,  rePONl  variant 

o 

G2E6  was  used  for  crystallization.  Crystals  of  rePONl -192Q  diffracting  to  2. 3 A  were  obtained 
as  for  G2E6  (192K)  (Ben-David  et  al,  2012;  Ben-David  et  al,  2013),  in  the  same  space  group  and 
with  very  similar  unit  cell  parameters  (Table  S2). 

Superimposition  of  rePONl- 192K  (PDB  code:  3SRE)  and  the  rePONl- 192Q  mutant 
structures  reveals  no  major  changes  in  backbone  or  side-chain  orientations  (Figure  S5A). 
However,  the  presumed  hydrogen-bonding  network  that  connects  H2  to  N 168  differs.  The  amide 

o  o 

group  of  Q192  is  in  close  contact  with  the  backbone  carbonyl  of  El 86  (2.9A,  versus  4.8A  to  the 

o 

e-amino  group  of  192K),  and  the  distance  to  the  side  chain  of  SI 66  becomes  shorter  (2. 6 A  vs. 
3.1A  in  192K;  Figure  S5).  F186  seems  to  be  involved  in  HDE  binding,  and  perturbations  of  the 
loop’s  backbone  at  this  residue  seem  to  have  a  marked  effect  as  indicated  by  the  F186G  mutation 
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(Table  SI).  It  thus  likely  that  the  K192Q  mutation  stabilizes  the  C-terminus  connecting  loop  in 
an  unfavorable  configuration,  such  that  the  HDL  binding  effects  on  H2  are  markedly  reduced. 

Interim  Conclusions 

The  mutational  and  structural  analyses  presented  here  provide  the  first  clear  evidence  for  the 
involvement  of  helix  H2  and  its  C-terminal  loop  in  the  binding  of  PONl  to  HDL,  and  in  the 
stabilization  and  stimulation  of  catalytic  activity  that  ensue.  The  residues  that  show  the  most 
pronounced  effects  appear  to  be  Y 185  and  FI  86,  which  reside  on  the  C-terminal  loop,  and  Y190 
and  W194,  which  reside  on  H2  itself  (Table  SI).  Thus,  our  refined  model  for  the  binding  of 
PONl  to  HDL  includes  three  sets  of  interactions:  the  N-terminal  membrane-anchoring  helix  HI, 
helix  H2,  and  H2’s  C-terminal  connecting  loop  (Figure  SI).  The  high  configurational  flexibility 
of  helices  HI  and  H2  is  likely  to  be  reduced  significantly  by  HDL  binding,  thus  accounting  for 
the  dramatic  stabilization  of  PONl’s  overall  configuration,  and  specifically  of  its  active  site  and 
catalytic  Ca"^^. 

At  the  technical  level,  the  cell- free  expression  of  PONl  in  the  presence  of  reconstituted  HDL 
particles  offers  a  new  approach  to  study  the  interactions  of  other  enzymes,  or  proteins,  with 
HDL,  or  with  lipoprotein  particles  in  general.  This  method  is  relatively  fast  and,  foremost,  allows 
production  of  stable  membrane-associated  enzymes  without  the  introduction  of  detergents. 

The  mechanism  of  catalytic  stimulation  of  PON  1  unraveled  here  seems  to  involve  a  complex 
network  of  interactions  that  spans  from  H2,  and  especially  its  C-terminal  loop  (residues  183- 
187),  into  the  active  site.  This  network  includes  two  key  active-site  residues,  D183  and  H184. 
These,  and  especially  D183,  are  entirely  conserved  (including  in  bacterial  PONs  that  are  >70% 
diverged  relative  to  mammalian  PONl).  Our  previous  studies  indicated  that  residues  183-184 
play  no  direct  role  in  catalysis,  but  why  these  positions  are  so  essential  had  remained  unclear.  It 
now  seems  that  they  play  a  primary  role  in  aligning  N168,  which,  in  turn,  ligates  the  catalytic 
Ca^^.  The  computational  simulations  support  this  model.  Specifically,  we  observed  that  changes 
in  H2’s  configuration  produce  changes  in  the  positioning  of  the  D183-H184  dyad  that,  in  turn, 
alter  the  interaction  of  N168  with  the  catalytic  Ca^^,  and  reduce  the  enzymatic  activation  energy. 
The  simulations  also  explain  a  critical  observation  regarding  the  selectivity  of  HDL  binding 
stimulation:  only  PONl’s  native  lactonase  activity  is  stimulated,  whilst  the  promiscuous 
paraoxonase  activity  remains  largely  unaffected. 
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Our  model,  and  the  newly  obtained  crystal  structure  of  rePONl-192Q,  also  account  for  the 
reduced  HDL  binding,  stability  and  stimulation  of  human  PON1-192Q  relative  to  the  192R 
isozyme.  The  alterations  to  the  hydrogen  bonding  network  include,  specifically,  the  appearance 
of  a  new  interaction  with  F186’s  backbone  carbonyl,  and  the  tightening  of  the  interaction  with 
S 166.  It  is  clear  that  the  structures  of  PONl,  whereby  the  enzyme  was  purified  and  crystallized 
in  detergent,  do  not  reflect  PONl’s  catalytically  optimal  configuration,  which  is  only  obtained 
upon  binding  to  HDL.  Thus,  although  no  changes  in  backbone  are  observed  with  192Q  and  R 
structures,  the  former  may  stabilize  a  suboptimal  configuration.  Indeed,  MD  simulations  of  the 
computational  models  of  both  the  192Q  and  192R  polymorphs,  along  with  rePONl’s  192K, 
indicated  that,  unlike  192Q,  192K/R  interact  extensively  with  both  D183  and  S166  (Sanan  et  al, 
2010). 

Finally,  our  results  also  rule  out  the  alternative  hypothesis,  that  Y71  mediates  the  effects  of 
HDL  binding.  This  despite  the  fact  that  Y71’s  hydroxyl  appears  to  be  H-bonded  to  D183  (Figure 
3),  and  the  cross-linking  data  that  implicated  it  in  HDL  binding  (Huang  et  al,  2013).  The 
strongest  point  that  contributes  to  elimination  of  this  hypothesis  is  that  the  Y71G  mutant  shows 
stability  comparable  to  that  of  wild  type  rePONl,  and  only  marginally  lower  enzymatic 
stimulation  (Table  SI).  Interestingly,  analysis  of  this  mutant  in  the  presence  of  a  detergent, 
tergitol,  showed  that  while  it  had  only  a  mild  effect  on  the  lactonase  activity  (~5-fold  decrease 
with  TBBL),  the  paraoxonase  activity  was  markedly  reduced  (~  120-fold  decrease).  Thus, 
although  the  highly  flexible  active-site  loop  at  the  tip  of  which  Y71  resides  plays  a  major  role  in 
substrate  specificity  (Ben-David  et  al,  2012),  its  effects  seem  to  be  HDL-independent. 

This  study  also  augments  our  limited  knowledge  of  how  lipid  interactions  modulate  enzyme 
structure,  catalytic  efficiency  and  specificity.  The  enzymes  best  studied  that  show  significant 
alteration  of  their  activity  upon  binding  to  lipid  membranes  (interfacial  activation)  are  lipases, 
and  especially  phospholipase  A2  (PLA2).  The  latter  shares  several  features  with  PONl,  such  as 
the  proximity  of  the  active  site  to  the  membrane  interface,  and  a  highly  mobile  active-site  loop. 
Analogously  to  the  findings  for  PONl,  truncation  of  the  N-  terminal  helix  of  PLA2  (Qin  et  al, 
2004),  or  even  point  mutations,  resulted  in  a  significant  decrease  in  both  PLA2’s  stability  and 
activity  (Chiou  et  al,  2008a;  Chiou  et  al,  2008b).  Also  in  resemblance  to  PONl,  the  N-terminal 
helix  of  PLA2  (residues  1-12)  shows  high  mobility  in  the  lipid-free  form  (van  den  Berg  et  al, 
1995).  This  configurational  flexibility  infiltrates  to  the  active  site,  and  specifically  to  the  catalytic 
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residues  H48  and  D99  (van  den  Berg  et  al,  1995).  It  was,  therefore,  postulated  that  the  lipid- 
bound  PLA2  adopts  a  fixed  conformation  that  reinforces  the  hydrogen  bonds  between  catalytic 
residues,  resulting  in  efficient  catalysis. 

The  revealed  mechanism  of  PONl’s  interfacial  activation  is  similar,  although  not  in  atomic 
details,  since  PLA2  and  PONl  differ  in  both  their  folds  and  their  active-sites,  but  in  the  overall 
effect  of  reducing  the  enzyme’s  conformational  ensemble,  and  aligning  the  key  catalytic  residues 
(N168,  and  the  catalytic  Ca^^),  thus  producing  lower  activation  energies  and  higher  catalytic 
rates. 

Beyond  these  specific  cases  of  interfacial  activation,  it  appears  that  active-site  floppiness, 
i.e.,  excess  of  mobility  of  active-site  residues,  is  a  challenge  to  enzyme  catalysis.  The 
coexistence  of  multiple  sub-states  is  counter-productive  with  respect  to  rate  acceleration,  namely 
catalysis  of  the  chemical  step,  as  is  seen  in  the  effect  of  H2  mobility  on  PONl’s  enzymatic 
activation  energy.  Widening  of  the  conformational  ensemble  may  also  account  for  the  gradual 
loss  of  rate  acceleration  with  increasing  temperature,  a  loss  that  seems  to  occur  in  all 

enzymes,  including  thermophilic  ones,  and  occurs  within  the  range  of  temperatures  in  which  the 
folded  state  is  fully  maintained.  Active-site  floppiness  also  seems  to  underline  designed 
enzymes,  where  the  conformation  of  key  active-site  residues  is  not  sufficiently  defined.  In  the 
course  of  the  evolutionary  optimization  of  computationally  designed  enzymes,  such  alternative 
configurations  are  minimized,  but  not  completely  eliminated  (Elias  et  al,  2014;  Giger  et  al,  2013; 
Khersonsky  et  al,  2012).  Indeed,  what  designed  active  sites  seem  to  be  lacking  most  is  the 
networks  of  interactions  that  align  catalytic  residues  such  N168  and  the  catalytic  Ca^^  in  PONl. 
These  networks  include  second-shell  residues  {e.g.  D183-H184)  and  third-shell  residues  (K192, 
SI 66  and  H2  residues;  Figure  S3).  The  essentiality  of  this  network  is  manifested  in  the  fact  that 
disturbance  of  even  its  most  far-reaching  periphery  (absence  of  HDL  interactions,  or  mutations 
in  192)  results  in  rate  declines  of  up  to  20-fold,  and  disturbance  of  second-shell  residues  such  as 
D183  leads  to  complete  loss  of  activity. 
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Table  SI:  Stimulation  of  activity  and  rates  of  inactivation  of  wild- type-like  rePONl  and  its 
mutants. 


Variant 

Activity  stimulation* 

Inactivation^ 

^app(pM)^ 

v^A%r 

V  IK 

^  max'  ^*^app 

(pM‘) 

Residual  lactonase 
activity  (%/ 

rePONl 

0.58±0.07 

(1) 

406±11 

(1) 

700±86 

(1) 

80±7 

(1) 

A20-rePONl 

>10* 

(>il7) 

— 

— 

20±6 

(i4) 

Y185A 

0.68±0.1 
(i  1.2) 

484±19 

(t  1.2) 

712±105 

(1) 

71±4 

(i  1.1) 

F186A 

2.4±0.3 

(i4) 

331±15 

(i  1.2) 

138±18 

(i5) 

11±1 

(i7) 

Y185A/F186A 

2.9±0.7 

(i5) 

183±17 

(i  2.2) 

63±16 

(i  11) 

6±2 
(i  13) 

Y185G 

>10* 

(>il7) 

— 

— 

15±3 

(i5) 

F186G 

1.5±0.4 

(i  2.6) 

159±13 

(i  2.6) 

106±30 
(i  6.6) 

7±1 

(i  11) 

Y190A 

1.3±0.3 

(i  2.2) 

268±17 
(i  1.5) 

206±49 
(i  3.4) 

59±2 
(i  1.4) 

W194A 

0.46±0.13 
(t  1.3) 

188±10 
(i  2.2) 

409±118 
(i  1.7) 

73±5 

(i  1.1) 

Y190A/W194A 

3±1 
(i  5.2) 

152±18 
(i  2.7) 

51±18 
(i  13.7) 

36±1 

(i  2.2) 

L198A 

2.2±0.2 
(i  3.8) 

518±19 
(t  1.3) 

259±23 
(i  2.7) 

47±6 
(i  1.7) 

Y71F 

0.8±0.1 
(i  1.4) 

467±14 

(t  1.2) 

584±75 

(i  1.2) 

75±5 

(i  1.1) 

Y71G 

1.70±0.09 
(i  2.9) 

430±7 

(t  1.1) 

253±14 

(i  2.8) 

73±6 

(i  1.1) 

H184Q 

1.5±0.2 

(i  2.6) 

367±12 

(i  1.1) 

245±34 
(i  2.9) 

34±3 
(i  2.4) 

H184T 

1.00±0.09 
(i  1.7) 

613±14 
(t  1.5) 

613±57 

(i  1.1) 

96±8 

(t  1.2) 

“  The  TBBL  concentrations  used  in  the  activity  stimulation  and  inactivation  are  0.25mM  and  0.5mM, 
respectively.  *  Stimulation  of  the  lactonase  activity  of  rePONl  by  different  concentrations  of  rHDL  (0- 
lOpM).  ^  The  apparent  affinity  for  rHDL  stimulation  is  the  concentration  of  rHDL  that  results  in  half- 
maximal  stimulation  of  the  activity.  values  are  presented  as  the  percentage  relative  to  the  expression 
of  rePONl  in  the  absence  of  rHDL.  Inactivation  of  rePONl  by  EDTA  (15mM).  ^The  residual  lactonase 
activity  is  presented  as  the  percentage  at  72  hr  of  incubation  at  37°C  relative  to  time  zero. 

*The  reaction  rates  did  not  show  saturation,  and  an  upper  limit  for  the  value  was  provided  based  on  the 
maximal  rHDL  concentration  used.  In  brackets  is  the  fold  decrease  or  increase  relative  to  the  wild-type-like 
rePONl. 
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Table  S2:  Summary  of  data  collection  and  refinement  statistics 


Data  Collection 

K192Q 

o 

Resolution  range  (A)'' 

40-2.3  (2.34-2.3) 

Space  group 

P432,2 

o 

Unit  Cell  Dimensions  (A) 

a=b 

93.5 

c 

144 

Number  of  reflections  measured 

410,343 

Number  of  unique  reflections" 

29,085  (1424) 

R  ^ 

^sym 

0.125  (0.684) 

Completeness  (%)" 

100  (100) 

Redundancy" 

14.1 (13.9) 

<I>  /  <s  (I)  >" 

27.7  (4.3) 

“  Values  in  parentheses  are  for  the  highest-resolution  shell. 
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Figure  SI:  The  proposed  model  for  anchoring  of  PONl  to  HDL  based  on  the  crystal 
structure  of  rePONl  (PDB  code:  3SRE).  This  model  suggests  the  interactions  of  PONl 
helices  HI  and  H2  with  the  HDL  surface.  Hydrophobic  residues  that  reside  on  helix  H2,  and 
on  H2’s  C-terminal  connecting  loop,  and  are  proposed  to  mediate  anchoring  to  HDL  are 
shown  in  orange. 
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Figure  S2:  Structural  properties  of  HI  and  H2.  The  analysis  was  based  on  the  wild- type-like 
rePONl  structure  (PDB  code:  3SRE). 

A.  Superposition  of  the  normal  modes,  as  predicted  by  elNemo,  indicates  the  high  mobility  of 
helices  HI  and  H2  (in  red). 

B.  The  two  edge  normal  mode  configurations  of  H2  and  its  C-terminus  connecting  loop, 
including  the  side  chains  of  H2’s  C-terminus  loop  residues  Y185  and  FI 86,  and  the  adjacent 
active  site  residues,  the  D83-H184  dyad,  the  active-site  mobile  loop  residue  Y71,  and  N168 
that  ligates  the  catalytic  calcium  (green  sphere).  The  inset  lists  the  distances  of  main-  and 
side-chains  between  these  two  extreme  configurations. 

C.  A  close  view  of  the  interactions  of  the  observed  part  of  a  detergent  molecule,  n-dodecyl- 
beta-D-maltoside  molecule,  DDM  (green;  shown  is  the  sugar  moiety,  the  electron  density  for 
the  lipophilic  chain  of  DDM  could  not  be  defined).  Shown  also  are  the  adjacent  residues  of 
Hi  (grey)  and  of  H2’s  N-terminal  connecting  loop  (magenta). 

D.  View  of  the  crystal  packing  interactions  within  different  PONl  molecules  (green  and 
cyan)  within  unit  cells  of  rePONl .  Shown  is  the  packing  surface  that  involves  interactions 
between  the  adjacent  H2  helices  and  their  C-terminal  connecting  loops. 
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□ 


Figure  S3:  Interactions  of  residues  D183  and  H184  with  first-shell  active-site  residues  (PDB 
code:  3SRG).  Depicted  are  putative  hydrogen  bonds  of  D183  and  HI 84  to  H2  residues 
(orange)  and  to  its  C-terminal  loop  (salmon),  as  well  to  the  catalytic  Ca^^-ligating  residue, 
N168  (grey).  Shown  is  the  lactone  analogue  2-hydroxyquinoline  (2HQ,  in  yellow)  and  Y71 
that  resides  in  an  active-site  mobile  loop  (cyan).  Our  model  suggests  that  HDL  binding 
results  in  the  fixation  of  a  well-defined  conformation  of  H2  and  of  its  C-terminal  loop. 
Immobilization  of  these  relatively  floppy  parts  of  the  enzyme  (Figure  2),  results  in  the 
alignment  of  the  D183-H184  dyad,  of  N168,  and  thereby  of  the  catalytic  Ca^^,  thus 
accounting  for  the  dramatic  increase  in  PONl’s  stability  and  lactonase  activity  when  bound 
to  HDL. 
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Figure  S4:  The  effect  of  H2  adjustment  on  the  active-site  architecture  and  on  substrate 
catalysis.  Superimposition  of  the  computational  models  with  helix  H2  positions  0  {color- 1),  2 
{color -2)  and  4  {color-3),  in  the  presence  of  TBBL  (A)  and  paraoxon  (B)  Shown  are 
remarkable  changes  in  the  positioning  of  D183,  H184,  N168  and  N224,  induced  by  the 
adjustment  of  helix  H2,  highlighting  structural  differences  between  the  models  produced  with 
TBBL  and  with  paraoxon. 
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Figure  S5:  The  crystal  structure  of  rePONl-192Q.  A.  Overlay  of  rePONl-192K  (PDB  code: 
3SRE;  blue)  with  the  newly  solved  structure  of  rePONl-192Q  (orange).  B.  The  network  of 
putative  hydrogen  bonds  spanning  from  H2  and  its  C-terminal  loop,  through  D183,  to  the 
catalytic  Ca^^-ligating  residue  N168  in  rePONl-192Q.  C.  The  same  analysis  for  rePONl- 
192K. 
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